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ABSTRACT 

In the non-linear phase of a dynamo process, the back-reaction of the mag- 
netic field upon the turbulent motion results in a decrease of the turbulence level 
and therefore in a suppression of both the magnetic field amplification (the a- 
quenching effect) and the turbulent magnetic diffusivity (the //-quenching effect). 
While the former has been widely explored, the effects of 77-quenching in the 
magnetic field evolution have rarely been considered. In this work we investi- 
gate the role of the suppression of diffusivity in a flux-transport solar dynamo 
model that also includes a non-linear a quenching term. Our results indicate 
that, although for a-quenching the dependence of the magnetic field amplifica- 
tion with the quenching factor is nearly linear, the magnetic field response to 
?7-quenching is non-linear and spatially non-uniform. We have found that the 
magnetic field can be locally amplified in this case, forming long-lived structures 



^The major part of this work was done during Gustavo Guerrero's visiting appointment at High Ahitude 
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whose maximum amplitude can be up to ~ 2.5 times larger at the tachocline 
and up to ~ 2 times larger at the center of the convection zone than in models 
without quenching. However, this amplification leads to unobservable effects and 
to a worse distribution of the magnetic field in the butterfiy diagram. Since the 
dynamo cycle period increases when the efficiency of the quenching increases, 
we have also explored whether the T^-quenching can cause a diffusion-dominated 
model to drift into an advection-dominated regime. We have found that models 
undergoing a large suppression in rj produce a strong segregation of magnetic 
fields that may lead to unsteady dynamo-oscillations. On the other hand, an 
initially diffusion-dominated model undergoing a small suppression in rj remains 
in the diffusion-dominated regime. 

Subject headings: MHD — sun: magnetic fields 



1. INTRODUCTION 



Over the past half a century, since the development of the first solar dynamo model by 
Parker (1955), significant investigations have been performed to find the saturation mech- 
anism that would limit the growth of a dynamo. Such mechanisms are likely to include 
feedback processes, such as the back-reaction of magnetic fields on the fiow fields, including 
mean fiows (differential rotation and meridional circulation), as well as turbulent flows. 

The feedback process that has been most extensively studied in the mean-fleld elec- 
trodynamics is the back-reaction of magnetic flelds on the helical part of the turbulent 
flow. This process, often known as a-quenching, was flrst believed to be mainly due to 
the back-reaction of the magnetic fleld on the convection, causing the suppression of ki- 
netic helicity, and hence quenching the inducing effects of the turbulent electromotive force 
Stix I (119721 ). Currently it is believed that the saturation process in the non-kinematic 
regime occurs due to the reduction of the kinetic a-effect caused by a magnetic contri- 
bution of opposite sign, coming from the equation of the a. This contribution appears 
asa product of the magnetic helicity conservation constraint. There exists a large liter- 



this topic: Kraichnan ( 


1979): 


Cattaneo k Vainshtein 


(1991^ 


: Gruzinov k Diamond 


1994 


); 

): 


Bhattachariee & Yuan 


(1995) 


: Cattaneo k Huehes Jl996h: 


Brandenburg k Donner 


1997 


Field. Blackman k Chou (1999): Blackman k Field 1 


2001) 


Field k Blackman (200 


2): Brandenburg k S 


(2005). 



However, the back-reaction due to induced magnetic fleld on the mirror-symmetric non- 
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helical part of the turbulent flow is a relatively less explored subject. Such back reaction has 
the effect of reducing the ed dy diffusivity. This back-r eaction process has been named the 
77-quenching, flrst derived by [Roberts fc Soward I (119751 ) using mean-fleld electrodynamics. 



Later, the role of ?7-quenching in a dynamo has 



( iRiidiger. Kitchatinov. Kiiker fc Schultz 



1994 



Tobias 



jeen investigated by several authors 



19961). primari l y in th e context of aQ- 



type stellar dynamo models. IRiidiger. Kitchatinov. Kiiker fc Schultz I (119941 ) incorporated 77- 
quenching in one- dimensional afl type and two-dimensional a'^fl type stella r dyna mo models. 
For supercritical dynamo regimes, IRiidiger. Kitchatinov. Kiiker fc Schultz I (119941 ) found that 
in one-dimensional dynamo models with r^-quenching, the fleld strength increased a bit, but 
not much, whereas the cycle period decreased signiflcantly. By contrast, for the ir two- 
dimensional stellar dynamo models, IRiidiger. Kitchatinov. Kiiker &: Schultz I (|l994l ) found 
a signiflcant fleld ampliflcation, almost two times more magnetic fleld was produced with 
?7-quenching, but the cycle period did not change much from that without r^-quenching. 
While fleld ampliflcation due to r^-quenching is intuitively expected, the change in dynamo 
cycle period, T, will depend on how T is determined in different classes of dynamo models. 
For example, in the convection zone aQ dynamo models, the cycle frequency {u) follows 
uj oc Q^-^i]^-^; hence an expected increase in cycle period with r^-quenching. 



By employing an aQ type interface dynamo model, iTobias I (119961 ) found that for weak 
magnetic flelds, the diffusivity is not much quenched, and the dynamo solutions are not 
i nfluen ced by the presence of r^-quenching. In the case of strong magnetic flelds, iTobias 



( 119961 ) showed that the diffusivity near the base of the convection zone can be so heavily 
quenched that the flelds can be trapped there without making their buoyant escape towards 
the solar surface. 



In a more recent calculation iGilman fc Rempel I (120051 ) showed, by solving the induction 
equation for the toroidal magnetic fleld component including the back-reactions of magnetic 
flelds on the turbulent diffusivity as well as on the shear, that there exists a competition 
between a fleld ampliflcation by r^-quenching and a fleld reduction due to the fact that the 
the Lorentz force fee dback on the shear moves the latitudinal shear layer away from the mid- 
latitudes. However, iGilman fc Rempel I (120051 ) argued that a signiflcant fleld ampliflcation 
might be possible if the latitudinal shear is replenished in a much shorter time-scale compared 
to the solar cycle. 

There is an effect common to all the above studies which either solved an aQ convection 



zone dynamo (IRiidiger. Kitchatinov. Kiiker fc Schultz Ill994l ). an interfa ce dynamo (iTobias 
1996h . or an induction equation for the toroidal magnetic fleld component (IGilman fc Rempel 
2OO5I ) - there is fleld ampliflcation due to r^-quenching. But the influence of r^-quenching on 
solar cycle features, namely the butterfly diagram and the evolutionary pattern of magnetic 
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fields in the convection zone and tachocline, has not yet been explored. 

Our aim here is to simulate a Babcock-Leighton flux-transport dynamo including the 
back-reactions of magnetic flelds on both the helical and non-helical parts of the turbulent 
flow, i.e. by including both the a-quenching and r^-quenching. We specifically seek the 
answers to the following questions: (i) what is a characteristic value of field amplification 
due to ?7-quenching in a Babcock-Leighton fiux-transport dynamo model? (ii) How is the 
butterfiy diagram changed or modified due to r^-quenching? (iii) Where in the solution 
domain that extends from pole-to-equator in latitude and from the tachocline to the solar 
surface in radial extent does 77-quenching have the most effect? (iv) Is the conveyor-belt 
mechanism preserved, or does it break down if rj is more and more quenched? (v) To what 
extent can the //-quenching affect the dynamo cycle period which is primarily determined by 
the meridional circulation in this class of models? (vi) Since the a-quenching saturates the 
growth of the dynamo field and the r^-quenching works in amplifying the field, how does the 
dynamo behave in presence of the competition between these two quenching mechanisms? 
(vii) Can the 77-quenching take a diffusion-dominated dynamo into the advect ion-dominated 
regime? 

In the next section we present the formulation of the model, including the prescriptions 
of the dynamo ingredients, such as the velocity fields (differential rotation and meridional 
circulation), diffusivity profile with r^-quenching formula, and the a-effect profile. We present 
the solution method, boundary conditions and initial conditions in section 3 and the detailed 
description of our results in section 4. We conclude in section 5. 



2. MATHEMATICAL FORMULATION 

In the mean field approximation, the MHD induction equation that governs the evolution 
of the large scale magnetic field is: 



= Vx[UxB + £ 

ot 



VtV X B], 



where B is the mean magnetic field and U is the large scale velocity field, rjT is the magnetic 
diffusivity and £ is the electromotive force, u x b, that represents the contribution of the 
small scale fluctuations upon the large scales. There have been several prev ious works in 



the literature where the effects of the latter term have been studied in detail flKiiker 



2OOII : iBonnano et al. 1 120021 : iKapyla et al. I l2006bl : I Guerrero fc de Gouveia Dal Pino 



et al. 



20081). 



These studies include the contribution in the mean-fleld dynamo equation explicitly from 
different parts of magnetic diffusivity tensor arising from small-scale turbulence along and 
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perpendicular to the rotation axis. In the present work since we are concerned in exploring 
the effects of the r^-quenching, we are going to neglect these detailed effects of £ other than 
the a-effect and the isotropic turbulent diffusivity. 

Working in spherical coordinates and assuming spherical symmetry, we can write Bp= V x 
{Ae^) and Be^, as the poloidal and toroidal components of the magnetic field, respectively, 
and considering Up and Or sin 9 as the meridional velocity field and the differential rotation, 
respectively, then we can split eq. ([T]) in the following pair of coupled partial differential 
equations for A and B: 



dA 
'dt 



+ -[up ■ V]isA) =nT{V^-^]A + Siir,9, B) 



dB 1 
dt r 



^{rUrB) + ^{ueB) 



- [Vr]T X (V X Be^)L + r/r V 



B 



(2) 
(3) 



where s = rsinO. The other terms will be described in detail in the next paragraphs. 



2.1. The velocity field 



The velocity field is one of the most important ingredients in the kinematic solar dy- 
namo. In the axisymmetric regime of a Babcock-Leighton flux-transport dynamo, the ve- 
locity field can be split into two large-scale mean-flow components in the azimuthal (0) 
and the meridional (r, 9) directions. The azimuthal component is the differential rotation, 
which is the responsible for the generation of the toroidal fields. The (r, 9) component is the 
meridional flow which transports the magnetic flux first poleward at the surface and then 
equatorward at the bottom of the convection zone, as in a conv eyor belt. This ingredient 



also plays a crucial role in determining the dynamo cycle period (jWang fc Sheeley Jr Ill991 



Dikpati fc Charbonneau 



1999 



KukeretaL 



200ll ) and the memory of the Sun's past mag- 



netic fields (jPikpati. de Toma fc Gilman II2006I ) if the advection dominates over the diffusion 
in determining the characteristic time-scale of the system. 

Recent observational developments gave us the access to an accurate a nd detailed profile 
of the differential rotation in the entire convection zone (see, for example, [Thompson et al. 
( 120031 ). and references therein). Observation of the meridional flow is a more difficult task. 
We now know its approximate shape near the solar surface. The temporal analysis of these 
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profiles support the use of tlie kinematic approximation for the solar dynamo, since the 
temporal variations in these profiles are small. From the observational point of view we can 
argue that the changes in these plasma movements due to the back reaction of magnetic 
fields are either small, or replenished in a time-scale much shorter than the solar cycle. 

In the pr esent calculation we inco r porat e the velocity profiles as prescribed in eqs. 
(4) and (5) of iDikpati fc Charbonneau I (119991 ). which resemble the helioseismology results 
regarding the differential rotation, and assume one meridional flow cell per meridional quad- 
rant, which is a likely assumption when both the observed all-latitude poleward flux and 
the mass conservation law are considered together. We notice that, despite several attempts 
have been made in fiux-transport dynamo m odels to include more than one convective cell 
(IBonnano et al. II2006I : iJouve &: Brun 1120071 ). no inferred magnetic field distribution agrees 
bette r with the observations than the one that is obtained by conside ring only one cell pat- 
tern (IDikpati et al. II2004I : iGuerrero fc de Gouveia Dal Pino Il2007bl ). Note however, that 
the depth of penetration, as well as the magnitude of the flow in the deeper layers are still 
uncertain in these models. In fact, nothing is known from observations regarding the struc- 
ture and the amplitude of the flow pattern inside the convection zone, except only very near 
the surface. 






-/ Ro 



atilude (degrees) 



Fig. 1. — Profiles of the main ingredients of the solar dynamo. The left panel shows the 
contours of iso-rotation (dotted lines) together with the meridional flow streamlines (solid 
contours), the dashed line at r=O.715i?0 indicates the center of the overshoot region. The 
middle panel shows the radial variation of the BL a term (solid line) and the magnetic 
diffusivity, riT (dotted line). The values for a are normalized to its maximum value ao, and 
the values of 77^ are at the right axis. The right panel shows the latitudinal profile of the BL 
a term. 

Both the contours of the iso-rotation and the streamlines of the meridional flow are 
shown in the left panel of Fig. 1. For all simulations presented below we consider a merid- 
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ional flow amplitude of uo/2=15 m and a depth of penetration rp=O.7i?0. Note that 
Uq/2 gives the maximum flow-speed (i.e. 15 ms~^) for a maximum uq = 30ms~^ for the 
mathematical prescription of meridional circulation we are using. In this calculation, we 
c onsider that the tachocline is centered at 0.7Rq with a thickness uj=0.03Rq (see e. g. 



Guerrero fc de Gouveia Dal Pino I (l2007al )). 



2.2. Diffusivity profile 



The magnetic diffusivity is another important ingredient in dynamo models, but we 
know very little about the amplitude and profile of this ingredient in the solar interior. 
At the surface, it should be of the same order as the supergranular diffusion, and at the 
radiative layer and beneath it should attain molecular values. However we do not have 
a good idea what value it should have in the bulk of the convection zone. Recently the 
turbulent diffusivity for large-scale magneti c fields has been estimated using the so-c alled 
test field m ethod for isotropic turbulence (ISur. Brandenburg fc Subramanian 1 12008| ) and 



convection (IKapyla et al. II2009I ). According to these studies, the turbulent diffusion is of 
the same order of magnitude as the first order smoothing estimate. This is in agreement 
with mixing length arguments for the kinematic regime in the range of currently accessible 
Reynolds number. 



While the calculation by iKapyla et al. I (120091 ) is a substantial advance and the best 
available today, given the limitation of computer powers, it is still far from solar-like condi- 
tions in terms of Reynolds number and the stratification. Inclusion of more realistic solar 
conditions as well as the j x B back-reaction that allows the formation of intense flux tubes 
may reduce the effective diffusivity in the bulk of convection zone. Furthermore, previous 
flux-transport dynamo studies have shown that the magnetic diffusivity is required to be 
one order of magnitude lower in the convection zone than at the surface in order for the 
dynamo to operate in the advection-dominated regime. So, in this work we use a diffu- 
sivity profile that has been used previou sly in several works (e.g 
Guerrero fc de Gouveia Dal Pino I (l2007al )). as follows: 



Dikpati et al. I fl2002[ ): 



, Vcz 
Vrz + — 



1 + erf 



r — Tr 



di 



2 



1 + erf 



r - rci 



(4) 



where r]rz=^0^ cm^ s 



, Vc 



=5 X 10^° cm^ s 



and ?7s=10^^ cm^ s ^ correspond to the values 



of the diffusivity at the radiative, convective and near-surface layers, respectively. The 
transition from the radiative to the convective layers is located at r.rc=O.715i?0 (the overshoot 
interface), with di=0.015RQ, and the transition from the turbulent convective zone to the 
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(sub-surface) supergranular diffusion layer is at rrc=0.96RQ, with (ii=O.O3i?0 (see the dotted 
hne in the middle panel of Fig. 1). 



2.3. Formulation of dynamo equations with //-quenching 

As mentioned earlier, the main goal of this work is to study the quenching of the 
turbulent diffusivity due to the presence of strong magnetic fields. For this aim we will 
assume that it will affect only the toroidal fields, and replace rjT in the equation for B (eq. 
E]) by 7]-. 



V 



Vt 



1 + {B/B, 



(5) 



which is the same algebraic form for the r^-quenching used by iGilman fc Rempel I (120051) . In 
the equation above, Bg is the value of the magnetic field at which rjT begins to be quenched. 
In principle, the poloidal fields also could contribute to the diffusivity quenching. However, 
here we focus only on the influence of toroidal fields in the saturation mechanism of the the 
turbulent diffusivity for two reasons: (1) in most of the a-Q solar dynamos, the dynamo- 
generated toroidal fields are about a thousand times stronger than the poloidal fields; (2) 
the amplitude of poloidal fields generated in a Babcock-Leighton flux-transport dynamo is 
of the order of a few hundred Gauss, much below the value of the lowest quenching field 
strength selected for the study in this paper. In a more realistic situation, the diffusivity is 
a tensor and its quenching is bound to have some effect due to the presence of poloidal fields 
also, and must be explored in future. 

Before we re-derive the dynamo equations with the inclusion of //-quenching, we briefly 
discuss the issue regarding the choice of r^-quenching formula. Analytical studies consid- 
ering 3D turbulence often yield a quenching formula proportional to \B\ rather than B'^ 
(IKitchatinov. Pipin fc Riidiger Ill994l : iRogachevskii fc Kleeorin 1120011 ). 



In turbulent forced MHD simulations, lYousef et al. I (120031 ) found that the suppression 



of the magnetic diffusivity follows the form: r7~r7j'o/(l + (^(^/Beq)"^), where the value of a 
depends on the geometry of the initial magnetic field (i.e. whether it is helical or not). Kapyla 
and colleagues (see http : arxiv. or gabsOSlO. 2298) have also recently performed numerical 
simulation where turbulent diffusivity is quenched. However, their results are not able to 
make a clear conclusion about the functional form of quenching to be proportional to \B\ or 
B^. In the present case, we will adopt a similar fo rmulation, but with Bq as a free parameter 



of the model. The results of lYousef et al. I (120031 ) indicate that in the solar convection zone 
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TjT is quenched approximately as in eq. ([5]) for fiducial choices of the value of Bq. 

We note that, with the quenching incorporated, t] is not only a function of r, but also 
depends on B{r,6,t), and hence r]{r,6,t) is a function of r, 6 and t. Due to the additional 
dependence of ri{r, 9, t) on 6 and t, the component of Vrj x (V x Be^) gives rise to two new 
terms in the induction equation. 

Thus the Equation (3) becomes: 



dB 1 

dt r 



g^lrUrB] + -[ueB] 



dr] 
dr 



dB 
dr 



B 

r 



+ 



1 dr] 

^de 



Bcote 



+r] V 



dB 

'de 



(6) 



2.4. Babcock-Leighton a effect 

The source of poloidal fields (eq. [2]) in a Babcock-Leighton dynamo is the decay of the 
tilted bipolar magnetic regions (BMR's), which form a net surface dipole moment that drift 
towards the poles and eventually cause the reversal of the polar fields. The place where 
the magnetic flux tubes, which are responsible for the formation of the BMR's, develop 
is uncertain, but due to various reasons we can assume that the flux tubes are formed 
at or below the base of the convection zone. There exists a vast literature on the topic 
of rising flux tube simulations that produce the tilt and emergence p attern at the solar 



surfac e which are in good agreement with observations (see, for example, iD'Silva fc Howard 



(Il993l ): iFan. Fisher fc McClymont I (119941 )). The combination of low magnetic diffusivity 
and helioseismically obtained differential rotation helps the amplification of toroidal fields 
there. 

Furthermore, the rising flux tube simulations indicate that their eruption latitude and 
the tilt acquired during their buoyant rise through the convection zone fit best with surface 
observations if flux tubes with magnitudes between 5 x 10^ and 10^ G are produced at the 
bottom of the convection zone. 0. 



^See, however, a discussion of altern ative possibilities in iDikpati et al. I ([2002); iBrandenburg I (j2005l ): 
Guerrero fc de Gouveia Dal Pino I (|2008l ) 
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Estimating the Babcock-Leighton a-effect by computing the buoyant eruption of mag- 
netic flux tubes followed by their decay is beyond the scope of this paper. So, in order to 
capture the properties of the BMR's described above, we simply use the following Babcock- 
Leighton a effect profile (jPikpati et al. Il2004j ): 



S,{r,e,B)=a{r,e)fQ{Br:)B,^ 



(7) 



This term is non-local in B (see below), Br^ being the radial average of B between r=O.7i?0 
and r=O.72i?0; this is a simple way to erupt toroidal magnetic flux from the bottom of the 
convection zone to the place where the a-effect is operating. We incorporate the following 
radial and latitudinal dependence of the a-effect: 



erf 



X sin 9 cos 9 



r — r2 
1 



erf 



^^3 



1 _^ e'^i('"/'^"^) 



where r2=O.95i?0, r3=i?0, (i2='i3=O.Oli?0, and 7i=30. The amplitude of the poloidal source 
is determined by ao, for which we assume a fixed value of 50 cm s~^ for the first set of 
simulations below (from now onwards we will explicitly quote the value of ao only when we 
are using a value different from this). As can be seen in Fig. 1, the a-effect is distributed 
mainly at the low latitudes peaking around 45°; it is also concentrated above O.95-R0 since it 
is expected that the main polo i dal field component is for med at the sunspot latitudes near 
the surface (IWang et al. Ill989l : IWang fc Sheeley Jr Ill99ll ). 



The quenching of the poloidal source term (the second term of eq. [7]) is given by: 



fQiB) 



1 + 



B. 



(9) 



This term has the same algebraic form as the a-quenching term used in turbulent mean 
field dynamo models. Its function here is to saturate the growing of the poloidal fields 
when the toroidal field at the base of the convection zone is around Bq=10'^ G. Thus the 
model will produce toroidal magnetic fields in the expected range, as explained above. The 
amount of poloidal field that will be produced from this toroidal field is determined by 
ao. We not e that this is not the unique way of capturing the physics of buoyantly erupted 
flux tubes. iNandy &: Choudhuri I (120011 ) have replaced this non-local quenching term by a 
different buoyancy mechanism term parameterized from simulations of flux tubes. 
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Although the equations ([5]) and (Q have the same functional form, they operate in 
opposite ways. We will discuss that in detail in §4. On the one hand, smaller values of Bg 
in eq. will produce lower values of rj which will, in turn, help to increase the amplitude 
of the toroidal component, B. On the other hand, lower values of Bq in eq. 1^ will limit 
the growth of B up to values around Bq. Both these terms are the sources of non-linearity 
in the model. 



3. SOLUTION METHOD, BOUNDARY AND INITIAL CONDITIONS 

We solve equations ([2]), and for A and B and r], with the coordinates r and 9 
covering the spatial range O.6i?0<r<i?0 and 0<6<7r/2, which spans from the outer most 
part of the radiative zone, through the overshoot layer (where the tachocline is located) and 
the entire convection zone, up to the surface, in the northern hemisphere. We have used a 
second order finite difference scheme for the spatial discretization; the Lax-Wendroff method 
for the first order derivatives, and centered finite difference for the second order derivatives. 



The te n aporal evolution is solved with th e ADI semi-implicit method (see lGuerrero &: MUnoz 



( 120041 ): iDikpati &: Charbonneau I (Il999l ). for details). In eq. ([5]), the dependence of rj with 
time is implicit, thus we update the value of 1] each half-time step with the previous values 
of B and then we calculate the derivatives of rj and use these values to solve the equation 
for B. 

The boundary conditions are: A = and i? = at the north pole {6 = 0); and at 
the equator B = 0, but A is coupled with the southern hemisphere in such a way as to 
ensure antisymmetric magnetic fields about the equator, so we demand ^ (r sin 9 A) = 0. 
At the bottom radial boundary we use A = -^(rB) = 0, and finally, at the upper radial 
boundary r = Rq, we consider a potential field boundary, i.e., B = and A coupled to 
an external vacuum field (V^ + 1/5^)^4 = 0. A complete de scription of the boundary con - 



ditions and their numerical imp lementation can be found in iDikpati fc Choudhuri I (jl994j ): 



Dikpati fc Charbonneau I (1199 



3). 



Since our intention is to explore how the 77-quenching affects a magnetic field that is well 
organized both in space and time, a plausible choice for an initial condition for the simulation 
is to start with a fully relaxed solution of the same dynamo without //-quenching. So we first 
obtain a fully converged solution by initializing the system with A = smO/r"^ if r > 0.715, 
y4 = otherwise, and 5 = 0, and allow it to evolve until 10^ yr with the r^-quenching in the 
Equation (6) turned off. 

Figure [2] shows the time- latitude, butterfly diagram for the last 60 years of evolution of 
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our reference simulation. It shows the gray scale (color) contours of the radially averaged 
toroidal magnetic field B^^ (in log-scale) with values above 12 kG, together with the contours 
of the radial field at the surface. It can be seen that the maximum amplitudes of the 
toroidal field are located below 45°, and satisfy the 7i/2 phase- lag observed with respect to 
the radial field. The maximum value of the toroidal magnetic field Br^ is ~ 22.9 kG, the 
maximum value of the radial field is 86.6 G, and the period of the entire cycle is 21.2 yr. 

Butterfly diagram 




0.0 15.0 50.0 45.0 60.0 

time (years) 

Fig. 2. — Butterfly diagram for the relaxed solution of the reference model. The dark 
(blue) contours represent positive toroidal fields and the light (red) contours represent the 
negative toroidal fields. The plotted values correspond to a radial average between O.7i?0 
and O.72i?0. Only contours for fields above 1.2 x 10^ G are plotted. The continuous and 
dashed lines represent the positive and negative radial fields at the surface, respectively. 
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0.0 0.2 0.4 0.6 0.8 1.0 

x/Ro 




0.0 0.2 0.4 0.6 0.8 1.0 

x/Ro 



Fig. 3. — Snapshots of the relaxed solution of our reference model at 0, T/8, T/4 and T/2, 
where T is the full period of the cycle. The gray scale (color) contours follow the same 
description as in Fig. [21 but in this case the line contours correspond to the total poloidal 
field. 
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The four frames of Fig [3] (a-d) show the same temporal evolution in the meridional cut. 
We see that the strongest toroidal fields begin to form at mid-latitudes inside the convection 
zone (r~0.8i?Q). The toroidal field penetrates slightly in the overshoot layer only at lower 
latitudes, but it is not substantially amplified there by the radial shear. This happens due 
to two reasons: first, since the meridional flow is not allowed in this model to go deep 
inside the overshoot layer, not enough poloidal fields (which are actually the source for the 
toroidal fields) can reach the radial shear layer; second, t he radial component of the poloidal 



field is much weaker there than the latitudinal component (IGuerrero fc de Gouveia Dal Pino 



2007al ). Fig [3] also reveals that the positive and negative poloidal magnetic fields (continuous 
and dashed lines, respectively) are produced at mid-latitudes at the surface and then migrate 
poleward following the plasma flow. 

We use the converged solution without the r^-quenching shown in Figures [2] and [3] as our 
initial {t = 0) configuration for our simulations with ?7-quenching and run the model for 200 
years more. In §4, we present our simulation results with 77-quenching. 
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4. Results 

4.1. The effect of a-quenching term in Babcock-Leighton dynamos 

Generally a-quenching is applied in most of the large-scale, mean-field dynamo models 
and the basic results are known. In a kinematic dynamo, the maximum value that the 
toroidal field can reach depends on the amount of poloidal field being generated by the 
Babcock-Leighton a-effect, and the later depends on the values of ao and Bq in Equations (IHl) 
and In order to explore the infiuence of a-quenching in more detail and to compare this 
infiuence with that obtained by implementing the ?7-quenching, we study how the maximum 
toroidal fields produced at different latitudes vary with the quenching field strength. The 
value of ao defines the non-dimensional number Ca=(yoRQ/Ticz=70 for a;o=50 cm s~^ and 
Vcz=5 X 10^° cm^ s~^, this value remains constant during all the simulations shown below. 
This guarantees that the dynamo efficiency CuCa is always the same. Then, we change the 
value of -Bo in eq. (Q between 5 x 10^ G and 5 x 10^ G and measure the maximum value 
that the toroidal field reaches at the numerical domain during one-half period. Intuitively 
we expect that the maximum value of B^ax should be larger for larger Bq. 

Fig. 4 presents B^ax as function of Bq for two different radii r=0.7RQ and r=O.8i?0, 
and three different latitudes, 10° (continuous line), 45° (dashed line) and 80° (dotted line). 
Figure 4 immediately reveals that maximum toroidal fields generated at different latitudes 
and different depths vary near linearly with the quenching field strength; the higher the 
quenching field strength the higher the dynamo-generated toroidal field. This means that 
the nonlinearity due to a-quenching is so weak even up to a quenching field strength of 10^ 
Gauss that the dynamo behaves virtually as if it is operating in the linear regime. 

We also see that for lower latitudes (< 45°), the maximum toroidal field is of the same 
order both in the convection zone (r=O.8-R0) and at the tachocline (r=O.7-R0). This is 
because the largest toroidal fields in a Babcock-Leighton fiux-transport dynamo are pro- 



duced mainly by the latitudinal shear 
than by the action of tachocline radia 
obtained previously by some authors ( 


working on the latitudinal poloidal fields rather 
shear on radial fields. This reinforces the result 


Jlempel 


2006: 


Dikpati, de Toma & Gilman 


2006; 


Guerrero & de Gouveia Dal Pino 


2007a 


). With no advective transport below the base of 



the convection zone, very little poloidal field diffuses down there, and due to the thinness of 
the tachocline, an even smaller radial component of those poloidal fields is available to be 
sheared there. However, the situation can be different if the magnetic fields can be trans- 
ported downwards by overshooting or magnetic pumpin g at the lower latitudes and then 



acqui re further amplification at the radial shear layer (iGuerrero fc de Gouveia Dal Pino 



20081). 
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10^ 10^ 10^ 

Bo (G) 

Fig. 4. — Maximum toroidal magnetic field as function of Bq at r=0.7RQ and r=O.8i?0 (see 
the a quenching function, eq. Q)) for three different latitudes, 10° (continuous line, cross 
symbol), 45° (dashed line, triangle symbol) and 80° (dotted line, X symbol). In the bottom 
panel, the continuous and dashed lines are almost coincident. In each panel, a -Bmax = Bq 
curve has been superimposed to show the linear proportionality of maximum toroidal fields 
produced by a certain a-quenching field strength. 
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4.2. The effects of ry-quenching in magnetic field evolution 

It has already been noted in the context of one-dimensional aQ dynamo, two-dimensional 
a'^Q dynamo and interface dynamo that, the quenching of the magnetic diffusivity due to 
the back-reaction of magnetic fields is a possible mechanism for further field amplification. 
In this subsection, we first explore in detail the evolution of magnetic fields in a diffusively- 
quenched Babcock-Leighton flux-transport dynamo, starting from the 2D reference-state 
solution described in §3. Subsequently we will also present the quantitative estimation of 
field amplification and the change in cycle period due to 77-quenching. 

In the case of the cc-quenching study, we ran our simulation for 5 x lO^G < -Bq < 5 x lO^G. 
We consider same range for Bq for the 77-quenching study: 5 x lO^G < Bq < 5 x lO^G, and 
keep all the other parameters the same as in the reference model described in §3. Using a 
fully converged solution of the reference model as the initial condition, we switched on the 
77-quenching and ran the simulation for 200 years. At this time, the system reaches a new 
steady state with cyclic variations oi A, B and rj. 

It is not surprising that, the smaller the Bq, the faster the 77^ is quenched, leading to a 
significant increase in toroidal field. Figures 5 and 6 show the temporal evolution of both, 
the toroidal and poloidal fields in the pole-to-equator meridional-cut for four successive times 
within half a cycle (left panels) for two different representative values of Bg (10"^ G and 10^ 
G, respectively). These are the solutions after 200 years' evolution. Right panel shows the 
change of 77 due to the quenching action, for two different latitudes, 10° (solid line) and 45° 
(dashed line). The non-quenched profile has been plotted also for comparison (red dotted 
line) . 

The common features in both cases (see Figures 5 and 6) are that the model exhibits 
a decrease in the diffusivity at the places where the toroidal field acquires considerable 
amplitude. This suppression can be as large as three orders of magnitude, as we see in 
Figure 6 for = 10^ G, but this suppression of 77 is not uniform - neither along the radial 
direction nor in latitude. Comparing the left and right panels in each of Figures 5 and 6, 
we see that the peaks and valleys are anti-correlated with the spatial distribution of the 
toroidal field amplitudes. Since we are presenting here the converged solution, this change 
in T] repeats in successive cycles. 
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Fig. 5. — Left panel: Snapshots for 4 different times within a half period of the relaxed 
model for 5^=10'^ G. Right panel: Turbulent diffusivity without quenching (red dotted line), 
for 10° (continuous black line) and 45° (dashed blue line). 




Fig. 6. — The same as in Fig. [5l but for 5^=10^ G. 
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In both cases, strong gradients in diffusivity in latitude as well as in depth occur, leading 
to the more efficient field amplification in those regions and hence, the formation of small 
regions of concentrated toroidal magnetic fields. Comparing Figures 5 and 6 with Figure 3, 
we can see the enhancement in magnetic fields even in places where a run without quenching 
does not exhibit strong toroidal fields. 

It is not difficult to understand why the diffusivity quenching leads to small regions of 
fiux concentrations in the computation domain. The average diffusion time of the fields at 
each portion of the domain where t] is strongly suppressed is larger than the neighbouring 
domains where the quenching is less effective. The lifetime of the toroidal fields is several 
years larger there than in the regions where the quenching is not as effective. Thus the 
former can undergo a prolonged amplification by the Vfl terms and can reach larger values 
than the neighbouring regions that have larger rj . 

An interesting feature we note in Figure 6 is the formation of small scale magnetic 
patterns at the overshoot tachocline (see in the left panels). This is associated with a 
variation in i] in finer spatial scale compared to what we see in Figure 5 as a function of 
depth. This is a consequence of the non-linear coupling between Bg and 77. The smaller the 
Bq, the narrower the rj profile, and the more segregation of the magnetic fields is produced!^. 
Due to large suppression in t] at the overshoot tachocline regions, the spot-producing toroidal 
fiux remains more frozen there, particularly in the case of Figure 6. In the meantime, two 
competing processes are going on - the prolonged shearing by the differential rotation, and 
the equatorward advection due to meridional flow. Since this advection must work against 
the diffusion of the flelds, we can think of it as "dragging" or "pulling" the fields along, at 
some net speed that is smaller than the meridional fiow there. If the equatorward advective 
drag partially wins at certain portions of these fields, those portions get torn out from the 
large-scale part of the fields, and locally reconnect. This happens predominantly on the 
equatorial side of the large-scale fields, and so more and more fragmentation of field takes 
place, leading to formation of small-scale structures. 



4.3. Influence of 77-quenching on field amplification 

In order to quantify how effective the ?7-quenching is in producing strong toroidal fields, 
we have plotted in Fig. 7 the maximum value of the dynamo-generated toroidal field as a 
function of Bq. As in Fig. 4, we show the results at three different latitudes: 10° (solid line). 



^We have tested the model calculation using a larger grid resolution, such as 200 x 200 grid points and 
the results remain unchanged 
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45° (dashed line) and 80° (dotted line), for two different radii, namely for 0.7i?o, the center 
of the tachocline, and O.SRq, the lower convection zone. 

At the tachocline (upper panel of Fig. 7), the curves show an interesting behaviour: 
at latitudes close to the poles (dotted hne), the toroidal field increases with decreasing Bq 
(i.e., with increasing quenching). The amplification factor could be up to ~ 2.5. We can 

understand this, because the poloidal fields carried down to the tachocline by the meridional 
flow, undergo prolonged shearing by the strong radial differential rotation there. However, 
we obtain an apparently counter-intuitive result that for lower latitudes (solid lines), there is 
a decrease of B with the increased quenching in r] (decrease in Bg). We see a little decrease 
of B at mid-latitudes also as Bq decreases, but the effect is not so pronounced, and appears 
only for very efficient quenching factors. For the fields at 10°, the factor of decrease could 
be as large as ~ 2.5. 

At the center of the convection zone (bottom panel of Fig. 7), the results are different; 
the magnetic fields for both low and mid-latitudes have the same increase, by a factor as large 
as ~ 2 with respect to the case with no quenching. With Bg — 5 x 10^ G the toroidal fields 
can reach values above lO^G. In the middle of the convection zone there is no radial shear. 
It is the latitudinal shear that works on the poloidal fields, so we see similar amplification 
for mid-latitude and low-latitude fields. If the flux tubes formed from these strong toroidal 
fields produced at the middle of the convection zone due to strong ?7-quenching, rise to the 
surface, their orientation may not agree with Joy's law. So we do not know whether such a 
strong ?7-quenching is working in reality, or some other processes are inhibiting the formation 
of flux tubes in the middle of the convection zone. 
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b) r = 0.8 R 




10^ 10^ 10^ 

Bq (G) 

Fig. 7. — Maximum toroidal magnetic field as function of Bq (see the eq. ([5]) for three 
different latitudes, 10° (continuous line, cross symbol), 45° (dashed line, triangle symbol) 
and 80° (dotted line, X symbol). At the bottom panel, only the results for 10° and 45° are 
shown, the maximum values of B for 80° are smaller than 10^ G. 



- 23 - 




Fig. 8. — Schematic diagram shows how, with the help of downward advective transport, the 
high-latitude toroidal fields (cross section of which is represented by the red loop, direction 
being perpendicular to the plane of the paper) can access the shear layer (within gray doted 
lines) is more efficiently amplified when there is strong ?7-quenching. This can be compared 
to the toroidal fields without r^-quenching (cross-section is represented by the green loop, 
direction being perpendicular to that loop), which diffuse more efficiently. By contrast, the 
upwelling flow at low-latitudes makes the toroidal flelds stay away from the shear layer when 
diffusivity is more suppressed (see the low-latitude red- loop). 
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The contrasting feature of the field amplification at high and low latitudes in the 
tachocline (see the dotted line and solid line in the upper panel of Figure 7) occurs, re- 
spectively, due to the difference in advcctive transport at high and low latitudes. If it would 
have been due to the suppression of the diffusivity and hence a prolonged shearing effect, we 
would have expected the increase in field amplitude in all latitudes. Clearly the directional 
transport plays some role here. We can understand this by using the schematic diagram 
shown in Figure 8. The advective transport being downward at high-latitudes, the poloidal 
fields there can reach the shear layer, for both cases when there is no 77-quenching and when 
there is strong 77-quenching. But in the latter case, the toroidal fields, originated from the 
sheared poloidal fields there, undergo less diffusion and stay at the tachocline shear layer for 
longer time than in the former case. The high-latitude toroidal fields (which are normal to 
the red loop in the figure) undergo further amplification due to local feeding by new toroidal 
field lines that are created by shear due to the strong tachocline differential rotation, com- 
pared to the case without 77-quenching where the toroidal field undergoes more diffusion (see 
the high-latitude green loop which represents a section normal to a bunch of toroidal field 
lines without r] — quenching). As a consequence, with the increase in the 77-quenching, there 
is a systematic increase in the high-latitude toroidal field amplitude. 

On the other hand, the upwelling fiow near the equator pushes the low-latitude poloidal 
fields upward, away from the shear layer, making it difficult to create toroidal fields at 
such latitudes. However a question arises here: the low-latitude poloidal fields are advected 
upward always, no matter whether the 77-quenching is present or not, but why do the low- 
latitude toroidal fields decrease with the increased 77-quenching (with decreasing Bg), instead 
of being independent of Bq7 Again we take the help of the schematic diagram in Figure 8 
to explain this feature. This happens due to the combination of decrease in diffusivity and 
upward transport. In spite of the fact that the meridional flow at low latitudes always takes 
the poloidal flclds away from the radial shear layer, there is a small amount of toroidal 
field being produced there and another amount produced at the convection zone due to 
latitudinal shearing. Given that the only mechanism that is able to transport the toroidal 
field downwards, at low latitudes, is the diffusive transport, in the case without 77-quenching 
(green loop), the toroidal field produced in the convection zone expands and reaches a por- 
tion of the tachocline. There, it encounters the existing amount of toroidal field and thus, 
increases its magnitude. However, in the case with 77-quenching, the toroidal field produced 
in the convective zone remains confined to a small region (red loop) and the toroidal field at 
the tachocline is not effectively increased. Hence the decrease in field amplification at low 
latitudes happens with increased 77-quenching. 
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4.4. Influence of r^-quenching in butterfly diagram and cycle period 

The increase in high-latitude toroidal fields and decrease in low latitude toroidal fields in 
the tachocline regions, with smaller Bq, will infiuence the butterfiy diagram accordingly. We 
recall again that the contours of the toroidal field that appear in the diagrams are computed 
from the radial average over the overshoot region. 

Figure 9 shows several butterfiy diagrams for different values of Bq. It can be seen 
how the butterfly wings, that are predominantly concentrated within latitudes < 45° when 
Bg is large (less quenching), move to higher and higher latitudes as Bq decreases (more 
quenching). This happens for the same reason that low-latitude fields at the tachocline 
decrease whereas the high-latitude fields increase with enhanced ?7-quenching. If a very large 
suppression in the rj occurs, the butterfiy diagram produced from this model will not be in 
accordance with observations. The ob vious question arises whether we can estimate how 



much //-quenching should be expected. iBrandenburg I (120011 ) showed in a direct numerical 



simulation of dynamos that the saturation level of the magnetic field is often close to the 
equipartition field strength given by, B^q = ^/JI^'pv? in MKS unit, or = ^Xjl pv? in CGS 
unit, where u is the rms value of the turbulent velocity, /zq is the conversion factor between 
CGS and MKS unit and p is the plasma density. With approximate values of the turbulent 
velocity of 5000 cms~^ at the base of the convection zone and density of 0.2 gm/cc, the 
equipartition magnetic field comes out to be approximately 10^ Gauss. Assuming that the 
back-reaction of the magnetic fields to quench the r/ will not start until the equpartition 
field strength is reached, the butterfiy diagram will not depart much from observations. By 
contrast, in the present calculations large departures from the observed butterfiy diagram 
occur for Bq < 500Gauss (see Fig. 9d), but according to the above arguments, 77-quenching 
should not occur for fields so far below equipartition. 

Note also that in the bottom panel of Fig. 9 the toroidal field shows some small-scale 
structures, as we saw in the toroidal field patterns in Figure 6. We repeat here that the pri- 
mary reason for the formation of these fragmented, small-scale structures is the competition 
between the two transport effects; the fields tend to remain more frozen due to the lowering 
in the diffusivity while the equatorward advective drag is pulling them, eventually causing 
their fragmentation. In order to check whether these are merely the numerical effects due to 
resolution problem, we have performed three experiments, namely doubling the resolution 
(experiment #1), reducing the advective speed to half of the value used in the present paper 
(experiment #2) and doubling the advective speed (experiment #3). The results (figures not 
included) indicate that the fragmented, small-scale structures do not go away with the in- 
creased resolution; nor do they go away when the meridional fiow-speed is reduced, but those 
structures are almost gone when the fiow-speed is doubled, because advection wins the com- 
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petition in that case. Perhaps a minor contribution into those small-scale structures comes 
from the numerical effect due to model diffusivity reaching the limit of grid- diffusion, but 
our experiment #3 confirms our physical explanation that the formation of those small-scale 
structures are the consequences of the two aforementioned competing processes. 

With smaller values of Bq, it is possible to reach larger magnetic fields; however the 
system does not reach a steady state. 

The changes in the butterfly patterns indicate that these diagrams remain consistent 
with observations up to a certain increase in the r^-quenching, but the influence of an en- 
hanced ?7-quenching is to make the butterfly- diagrams depart further from the observations. 
The shift in the butterfly wings towards higher latitudes could cause another problem in the 
dynamo model, namely loss in the coupling between the two hemispheres across the equator 
and hence, a shift to the quadrupole parity in the solution when solved in a full spherical 
shell. In that case, to restore the observed dipolar parity of the large-scale solar magnetic 
fields, the help of additional downward transport using turbulent pumping may be required 
(Guerrero & de Gouveia Dal Pino 2008). We leave those studies for the future. 

In Fig. 9 we also see an increase in the dynamo cycle period with increasing ?7-quenching, 
i.e. with decreasing Bq. The number of butterfly wings produced in 60 yr decreases. Figure 
10 shows a plot of dynamo cycle period as function of Bq, and the decrease in cycle period 
with increased 77-quenching is very clear. This happens due to the competition between the 
flux-freezing effect by the smaller diffusion and advective drag due to meridional flow. The 
suppression in 77 due to quenching works against the advective transport of flux. 
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Fig. 9. — Butterfly diagrams for different values of Bg. Top: Bg=10^ G; second from top: 
Bg=10^ G; third from top: Bg=10^ and bottom: Bg=5 x 10^ G. 
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4.5. Advection-dominated versus diffusion-dominated dynamos 

In Babcock-Leighton dynamo models, the meridional flow is the conveyor belt that 
carries the magnetic field, both at the surface in order to produce a new poloidal field, and 
at the bottom of the convection zone where it transports the toroidal fields in the direction 
of the equator in such a way that this flow dominate s over other parameters in setting 



the period of the cycle (jPikpati fc Charbonneau Ill999l ). But for this process to occur the 



advective term in eq. [T] must dominate the diffusive term in determining the time-scale of 
the system. Thus this class of models require values of rjcz^'i x 10^^ cm^ s^^. However, 
the values inferred from the mixing length theory at the surface are one to two orders 
of magnitude larger than the values considered in the bulk of the convection zone in the 
advection-dominated flux-transport dynamo models. To our knowledge, until now there is 
no accurate estimation of the magnetic diffusivity at the convection zone and beneath. This 
limitation ha ve recently led t o crit ic isms to the flux-tran sport scenario (see more about this 



discussion in ICharbonneau I (120071 ): lYousef et al. I (120031 )). 



In previous sections, we have described the action of the r^-quenching in a Babcock- 
Leighton dynamo and found that, depending on the quenching parameter Bg, the diffusivity 
can be locally suppressed by up to three orders of magnitude and this effect can also increase 
the cycle period (Figure 10). This is an indication that the average radial value of rjcz is 
also being suppressed. One question arises here: will diffusion-dominated dynamos, which 
in general produce much faster cycles, change to advection-dominated dynamos due to the 
suppression of rj by quenching mechanism, and produce a cycle period similar to the observed 
sunspot cycle? 

We have performed simulations to try to answer this question. We have increased pro- 
gressively the diffusivity from 5 x 10^° cm^ s^^ (the value employed in the previous calcula- 
tions) to the largest allowed value, namely 5 x 10^^ cm^ s~^, for which we obtain an oscillatory 
solution (not decaying due to large diffusivity). The latter case is the diffusion-dominated 
regime, and the cycle period is small, determined by the diffusivity values. We then per- 
formed the same numerical experiments as before by switching on the ?7-quenching, for the 
two representative values of Bq, 10^ and 10^ G and maintained the constant dynamo efficiency 
deflned by CnCa, for all the simulations. We note that keeping CnCa = aoVLeqFl? /vIz = cte 
implies a change in as r/cz changes (see the Table 1 for the values used for for each 
simulation, as well as the maximum values for the average toroidal fleld at the base of the 
convection zone, B^^., and the radial fleld, B^., at the surface). 

In Figure 11, we present the dynamo cycle period as function of magnetic diffusivity 
for the two cases, with Bq = 10"^ G (pluses) and Bq = 10^ G (triangles). The models with 
Bq=10'^ G (plus symbols) present two different regimes for the slope of the curve of the 
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Fig. 10. — Half-cycle period as function of Bq. 

period versus rjcz- For values of ?7c^< 10^^ cm^ s^^, the cycle period does not vary much with 
the value of rj; it is primarily determined by the meridional flow speed. For ricz>W^^ cm^ 
s~^, the cycle period varies more rapidly with rjcz, indicating that the models are operating 
in the the diffusion-dominated regime. 

For the case of 5^=10'^ G, the highest value of rjcz, that allows a steady state solution 
with a well deflned period is 10^^ cm^ s~^. In this case, the diffusivity is highly intermittent 
with larger gradients than in all the previous calculations. However, we flnd from the plot 
of triangles in Figure 10, a smooth change of the period as function of rjcz- 

Figure 12 shows a butterfly diagram for one of these diffusion-dominated cases. For 
these models (with higher r]cz), we also see in Figure 12 that the phase difference between 
the toroidal fleld at the lower latitudes and the radial flelds near the poles is ~ tt rather than 



These results clearly reveal that a diffusion-dominated dynamo with meridional circu- 
lation remains in a diffusion dominated regime even if the diffusivity is locally suppressed 
by the back-reaction of magnetic flelds. However, we emphasize the fact that the amplitude 
and profile of the magnetic diffusivity inside the convection zone are still unknown. 



7r/2. 



- 30 - 



o 



1 r 



(D 8 - + ^ 

Q= : : 

ji) 6 - + ^ 

o : + = 

o 4 r T 

? 2 r T 

F ^ ^ ^ ^ ^ ^ : 

10^^ 10^^ 10^° 

7]^^ (cm^ s"'') 



Fig. 11. — Half-cycle period as function of rjcz for 5^=10'^ G (triangle symbol) and 5^=10^ 
G (cross symbol). 
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Table 1: Simulation parameters [Bg, rjcz and ao) and model-output (maximum toroidal field 
at convection zone base, maximum surface radial field and dynamo cycle period). 



BAG) 


r]cz (cm^ s"^) 


cto (cm s ^) 




Br^a^ (G) 


T/2 (yr) 


10^ 


5 X 10^° 


50 


4.8 X 10^ 


153.7 


13.65 


10^ 


6 X 10^° 


72 


5.0 X 10^ 


170.0 


13.23 


10^ 


7 X 10^° 


98 


7.5 X 10^ 


543.1 


12.75 


10^ 


8 X 10^° 


128 


7.9 X 10^ 


832.8 


12.44 


10^ 


1 X 10^^ 


200 


9.3 X 10^ 


981.428 


11.55 


10^ 


5 X 10^° 


50 


3.6 X 10^ 


121.9 


11.74 


10^ 


6 X 10^° 


72 


4.5 X 10^ 


174.1 


11.57 


10^ 


7 X 10^° 


98 


5.3 X 10^ 


252.2 


11.37 


10^ 


8 X 10^° 


128 


6.1 X 10^ 


345.2 


11.14 


10^ 


9 X 10^° 


162 


6.7 X 10^ 


452.0 


10.90 


10^ 


1 X 10^^ 


200 


7.2 X 10^ 


572.4 


10.64 


10^ 


2 X 10^^ 


800 


9.5 X 10^ 


2547.5 


8.07 


10^ 


3 X 10^^ 


1800 


1.0 X 10^ 


5998.7 


6.23 


10^ 


4 X 10^^ 


3200 


1.1 X 10^ 


10635.7 


5.07 


10^ 


5 X 10^^ 


5000 


1.2 X 10^ 


16515.5 


4.32 
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5. SUMMARY AND COMMENTS 

We have explored here the effects of diffusivity quenching on Babcock-Leighton flux- 
transport solar dynamo models. We used as initial condition a converged solution of a 
dynamo model that reproduces most of the main features of an observed solar butterfly 
diagram. The //-quenching was then included in the model through an algebraic function 
that is similar to the usual a-quenching formula. After some years of evolution, the system 
reaches a new steady state configuration in which the poloidal and toroidal components of 
the magnetic field, A and B, respectively, as well as the magnetic diffusivity, t], exhibit a 
cyclic behavior. 

With the new rj profile (Equation 5), the decrease in the diffusivity can be as large 
as three orders of magnitude, but it is not homogeneous over the whole domain since it 
presents a pattern with peaks and valleys that anti-correlate with the spatial distribution of 
the amplitude of the toroidal fields - the smaller the Bq (the value of the magnetic field at 
which the diffusivity begins to be quenched; eq. 5) the narrower the final rj profile (Figs. 
|6] and [5]). This spatial fluctuation in the magnetic diffusion results in the formation of 
small and long-lived regions of concentrated magnetic field that appear predominantly at 
the equatorial part of the overshoot tachocline as well as in the middle of the convection 
zone. The role of the meridional flow is very important in this result, since the toroidal 
fields have time enough to increase and be dragged along with the poloidal field lines before 
being dissipated. Note that if these strong flux tubes produced at the middle convection 
zone emerge to the surface they may not agree with the Joy's law. 

We have found that, contrary to the effects of a-quenching, which shows an almost linear 
coupling between the saturation fleld, Bq and the flnal value of a (see Fig. H]), the dependence 
between the diffusivity saturation fleld, Bq, and the flnal value of rj is predominantly non- 
linear, especially at the base of the convection zone where both the radial and the latitudinal 
shear are competing with the advective transport and with diffusive spreading. It was found 
that at the tachocline (r=O.7i?0) the magnetic fleld can be amplifled by a factor up to ~ 2.5 
at the highest latitudes, while for the equatorial regions the magnetic fleld decreases by 
approximately the same factor. On the other hand, at the center of the convection zone 
(r=O.8i?0), the magnetic fleld for both low and mid-latitudes have the same increase, by a 
factor as large as ~ 2 with respect to the no-quenching case. 

The consequence of these effects on the butterfly diagram is the increase (decrease) 
of the toroidal flelds at the high (low) latitudes and the increase of the cycle period for an 
enhanced 77-quenching. This new distribution of magnetic flelds with latitude not only shows 
a gradual departure from observations, but also could cause other problems in the dynamo 
models, such as the loss of coupling between hemispheres when solved in a full spherical 
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shell, resulting in a quadrupole parity solution. 

Since the period of the cycle diminishes when the efficiency of the quenching increases, 
we have also explored if the //-quenching can make a diffusion dominated dynamo, which is 
characterized by a small cycle-period, to evolve to an advection dominated dynamo. We have 
found that a diffusion-dominated dynamo with a meridional flow remains in the diffusion 
regime even when the diffusivity is locally suppressed due to the strong magnetic fields. 

The results summarized above indicate that in the scenario of a pure Babcock-Leighton 
dynamo, with a meridional flow operating as a conveyor-belt and strong magnetic flux tubes 
emerging from the base of the convection zone, the turbulent diffusivity is probably weakly 
suppressed, i.e., it is quenched only for high values of the magnetic field. This implies 
that its role in the amplification of the magnetic field to values above the equipartition 
field should not be significant enough. Notice, however, that other effects, like turbu- 
l ent pumping, which have been demons trated to be important in the dynamo operation 



( iGuerrero fc de Gouveia Dal Pino II2008I ). have not been considered here. The contribution 
of turbulent pumping might, for instance, result in a different transport of the magnetic 
fields, changing the parts of the parameter space where the r^-quenching would become dom- 
inant. This will be explored in forthcoming work. Finally, it is also important to remark 
that the quenching effects of the diffusivity still need to be explored in other classes of dy- 
namo models, for example, the ones operating with a surface shear layer. These will also be 
considered in future work. 
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